The analyticity region of the hard sphere gas. Improved bounds 



Roberto Fernandez 

Laboratoire de Mathematiques Raphael Salem 
UMR 6085 CNRS-Universite de Rouen 
Avenue de l'Universite, BP.12 
F76801 Saint-Etienne-du-Rouvray (France) 

Aldo Procacci 

Departamento de Matematica UFMG 
30161-970 - Belo Horizonte - MG Brazil 

Benedetto Scoppola 

Dipartimento di Matematica Universita "Tor Vergata" di Roma 
V.le della ricerca scientifica, 00100 Roma, Italy 

February 1, 2008 



Abstract 

We find an improved estimate of the radius of analyticity of the pressure of the hard-sphere gas 
in d dimensions. The estimates are determined by the volume of multidimensional regions that can 
be numerically computed. For d — 2, for instance, our estimate is about 40% larger than the classical 
one. 

In a recent paper [3] two of us have shown that it is possible to improve the radius of convergence of 
the cluster expansion using a tree graph identity due to Penrose [2] , see also [SJ Section 3] . In this short 
letter we use the same idea to improve the estimates of the radius of analyticity of the pressure of the 
hard-sphere gas. 

The grand partition function S(z, A) of a gas of hard spheres of diameter R enclosed in a volume 
A C K is given by 

A ) = ^2 ~7 / dxi — dXn ex p{ _ U(xi-x 3 )} 

n=0 ' l<i<j<n 

with 



U(x -y) = 



i£\x-y\>R 
+oo if \x — y\ < R 



where \x — y\ denotes the euclidean distance between the sphere centers x and y. The corresponding 
pressure is limA^ R d P(z, A) (limit in van Hove sense), where 

P(z,A) = -LogS(z,A) 
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(|A| denotes the volume of the region A). The cluster expansion, in this setting, amounts to writing the 
preceding logarithm as the power series (see e.g. pQ) 



logSO, A) = ^2 — / dx 1 ...dx n ^2 (-1) 



\tl\ 



(1) 



where the graph g{x\, ....£„) has vertex set {1, ...,n} and edge set E(xi, ....:£„) = {{i, j} : \xi — Xj\ < i?} 
(that is, if the spheres centered at Xi and Xj intersect), G n is the set of all the connected graphs with 
vertex set {1, n}, and \g\ denotes the cardinality of the edge set of the graph g. Only families (x\, ....£„) 
for which g(x\, ....x n ) is connected contribute to (1); such families represent "clusters" of spheres. 

The standard way to estimate the radius of analyticity of the pressure is to obtain a A-independent 
lower bound of the radius of convergence of the series 



\P\ (z,A) 



|A| 4-| nl 



dx\...dx r , 



E (-D 



gCg(xi ,....x n ) 



(2) 



This strategy leads to the classical estimation (see e.g. j6], Section 4) that the pressure is analytic if 

1 



* < 



eV d (R) ' 

where Vd(R) is the volume of the d-dimcnsional sphere of radius R (excluded volume). 



(3) 



Our approach is based on a well known tree identity. Let us denote by T n the subset of G n formed 
by all tree graphs with vertex set {1, n}. Given a tree r G T n and a vertex i of r, we denote by dj 
the degree of the vertex i in r, i.e. the number of edges of r containing i. We regard the trees r G T n 
as rooted in the vertex 1. This determines the usual partial order of vertices in r by generations: If 
it, v are vertices of r, we write u -< v — and say that u precedes v — if the (unique) path from the root 
to v contains u. If {u, v} is an edge of r, then either v -< u or u -< v. Let {u, v} be an edge of r and 
assume without loss of generality that u -< v, then u is the called the predecessor and v the descendant. 
Every vertex v G r has a unique predecessor and s v = d v — 1 descendants, except the root that has no 
predecessor and s v = d v descendants. For each vertex v of r we denote by v' the unique predecessor of 
v and by v , . . . v s ' v the s v descendants of v. The number s v is called the branching factor; vertices with 
s v = are called end-points or "leaves" . 

Penrose [4] showed that the sum in (1) is equal, up to a sign, to a sum over trees satisfying certain 
constraints. We shall keep only the "single-vertex" constraints: descendants of a given sphere must be 
mutually non-intersecting. This implies that 



]T (-1)1*1 < J2 «v(*l>""Zn) 
tET„ 



(4) 



where 



{1 if \x v — x v / \ < R and \x v % — x v j | > R, Vv vertex of r, {i, j} C {1, . . . , s v }, 
( 5 ) 
otherwise 



Hence from (2) and (4) we get 



with 



|P|(*,A) < K^aM 



n=l t£T„ 



S ^( T ) — TTT / ^(^1, ....X n )dxi...dx T , 

l A l JA" 



(6) 
(7) 



2 



By (5) we have 

n 

Sa(t) < ff rf (di)n^R-l) (8) 

1=2 

where di is the degree of the vertex i in r, 

9d{k) = dxi...dx k = R dk / dyi...dy k 

J Kl<fl J i» 4 i<i 

|a=i-a=jl>« l»i-Vj!>l 

for fc positive integer, and <?d(0) = 1 by definition. It is convenient to write 

9d {k) = [V d (R)] k g d (k) (9) 

with 

Uk) = Wfiwf™* dyi '" dyk (10) 

\Vi— Vj |>1 

for & positive integer and §d(0) = 1. We observe that < 1 for all values of k. From (8)-(10) we 

conclude that 

n 

Sa(t) < [VdiR^g^Umm^Udi-l) 

n 

= [V d {R)) n - x g d {dx)Jlg d {di-l) ■ 

i=2 

The last identity follows from the fact that for every tree of n vertices, <£]. + ••• + d n = 2n — 2. The 
r-dependence of this last bound is only through the degree of the vertices, hence it leads, upon insertion 
in (6), to the inequality 



dl + ... + e£ Tl = 27i-2 



The last quotient of factorials is, precisely, the number of trees with n vertices and fixed degrees d%, . . . ,d n , 
according to Cayley formula. 

At this point we can bound the last sum by a power in an obvious manner. The convergence condition 
so obtained would already be an improvement over the classical estimate (3). We can, however, get an 
even better result through a trick used by two of us in [SJ. We multiply and divide by a™ -1 where a > 
is a parameter to be chosen in an optimal way. This leads us to the inequality 

\P\( ^ < a f- (M^W)" Mdi)a dl A &(<*«-!) a*" 1 

dx+...+d n =2n— 2 



< 



The last series converges if 



71=1 V ' \ S>0 



gd(s)a s 



\z\V d (R) < " 



C d {a) 
where 

9d(s) 



C d {a) = 



s! 

s>0 
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(this is a finite sum!). The pressure is, therefore, analytic if 

\z\V d (R) < max-^- . (11) 
a>0 Ld(a) 

This is our new condition. 

Let us show that for d = 2 the quantitative improvement given by this condition can be substantial. 
In this case 

Ca(o) = 2^—r a 

s=0 S - 

where, by definition, 32(0) =32(1) = 1- The factor 172(2) can be explicitly evaluated in terms of straight- 
forward integrals and we get 

52(2) = IT 

The other terms of the sum can be numerically evaluated using a simple Montecarlo simulation, obtaining 

52(3) = 0,0589 <? 2 (4) = 0,0013 52 (5) < 0,0001 

1 / 2 

(a value for which (j^hs is close to its maximum) we get 
\z\V 2 {R) < 0.5107. 



Choosing a 



3V3 



This should be compared with the bound \z\ V^(i?) < 0.36787 . . . obtained through the classical condition 
(3). 



Acknowledgments: We thank Sokol Ndreca for useful discussion concerning the numerical evaluation of 
the g's. A. P. thanks CNPq for financial support. 



References 

[1] D. Brydges: A Short Course on Cluster Expansion, Les Houches 1984, K. Osterwalder, R. Stora 
eds., North Holland Press (1986) 

[2] O. Penrose (1967): Convergence of fugacity expansions for classical systems. In "Statistical mechan- 
ics: foundations and applications", A. Bak (ed.), Benjamin, New York. 

[3] A. Procacci and B. Scoppola: Polymer gas approach to N-body lattice systems. J. Statist. Phys. 96, 
49-68 (1999) 

[4] R. Fernandez and A. Procacci: Cluster expansion for abstract polymer models. New bounds from an 
old approach, Preprint arXiv math-ph/0605041, to appear in Comm. Math. Phys. 

[5] Ch.E. Pfister: Large deviation and phase separation in the two-dimensional Ising model, Helv. Phys. 
Acta 64, 953-1054 (1991) 

[6] D. Ruclle (1969): Statistical mechanics: Rigorous results. W. A. Benjamin, Inc., New York- 
Amsterdam. 



4 



